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Abstract: We consider the Higgs Lepton Flavor Violating process h —> r/r, in which 
CMS found a 2.5 a excess of events, from a model independent perspective, and find that 
it is difficult to generate this operator without also obtaining a sizeable Wilson coefficient 
for the dipole operators responsible for tau radiative decay, constrained by BABAR to 
BR(r — > n'y ) < 4.4 x 10~ 8 . We then survey a set of representative models for new physics, 
to determine which ones are capable of evading this problem. We conclude that, should this 
measurement persist as a signal, type-III Two Higgs Doublet Models and Higgs portal-like 
models are favored, while SUSY and Composite Higgs models are unlikely to explain it. 
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1 Introduction 

Now that we know of the existence of a Higgs boson with a mass of 125.7 ± 0.4 GeV [1] 
the focus has shifted to the precise determination of its other properties, with the aim of 
determining if they are all consistent with the Standard Model (SM) Higgs boson. All the 
data collected so far is consistent with the predictions of the SM [2], but the uncertainties 
are large and new physics could still be hidden, specially in the couplings of the Higgs boson. 
Deviations on these couplings could point to physics Beyond the SM (BSM), and since the 
Higgs couples directly to most of the other SM fields, it is an excellent probe for a multitude 
of different BSM models. New fermions not far above the electroweak scale, for instance, 
could alter the effective couplings of the Higgs boson with gluons and photons, and other 
sources of electroweak symmetry breaking could move its couplings to SM fermions away 
from the predicted values. 

A lot of experimental effort has been put into measuring the decay rates both in the 
dominant and in subdominant channels (that are either forbidden at tree level or suppressed 
by small couplings). The subdominant decays are specially interesting, because the decay 
rates predicted by the SM are so small that they can be easily dominated by BSM physics, if 
any is present. In this paper, we will concentrate on Lepton Flavor Violating (LFV) decays 
of the Higgs: h —> Ulj (which is understood as the sum of h —> iflj and h —> ). These 

decays are indirectly constrained by low energy data such as muon and tau rare decays, the 
g — 2 and Eletric Dipole Moments (EDM) of electrons and muons. 

The rare decays of the muon are so well tested that the ensuing constraints put the 
LFV decays of the Higgs to eg, beyond the reach of the LHC [3]. The constraints on tau 
decays, on the other side, are more allowing and recent papers suggested that FV decays 
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of the Higgs into tau could be observed in the LHC [3-5] with BR(/i —> t^l) < 10 -1 . That 
possibility motivated a dedicated search for h —> t/j, events, the results of which have been 
recently released by CMS [6], with an excess of 2.5 a 1 . The conclusions of [3-5] were 
obtained in fairly model independent approaches, using effective Lagrangians and a small 
variety of operators. On the other side, many papers have explored Higgs FV in specific 
UV completions of the SM without any additional light particles [7-19] and some actually 
obtain, given the indirect constraints, a much smaller BR(/i —> t/j), below the current 
reach of the LHC [7, 10, 10-13, 18, 19]. It is this apparent contradiction between model 
independent approaches and specific UV continuations that motivates this paper. 

We will approach the problem, in section 2, from a model independent perspective, 
using an effective theory and including dimension six and eight dipole operators, which 
have been mostly neglected in previous works on LFV Higgs decays. We will argue that the 
inclusion of those operators changes how the bounds from low energy data can be translated 
to constraints on LFV Higgs couplings, making the constraints in general more restrictive. 
We will also show that, in many realistic and simple UV continuations, dipole operators 
are generated with sizeable Wilson coefficients, and it is natural that in those cases the 
h —> Tfi branching ratio will be much smaller than predicted by analyses that neglected 
those operators. 

In section 3 we will explore some representative classes of models to determine which 
ones generate which operators, and what are the relations between the Wilson coefficients. 
That classification will allow us to have a general view of the problem, separating models 
that predict branching ratios much smaller than the current experimental reach from those 
that are actually favored in case the 2.5 a excess, observed by CMS in the h —» r/i channel, 
gets more significant with increasing luminosity. 

2 Model Independent Considerations 
2.1 Effective Field Theory 

In the Effective Field Theory (EFT) approach, we consider the most general Lagrangian 
compatible with the SM local symmetries and containing no new degrees of freedom below 
a set scale A: 

£eft = 4m + T. , C 2 - 1 ) 

i> 4 

where 4 m is the renormalizable SM Lagrangian, and contains the higher dimensional 
operators (of dimension i) generated by integrating out heavy fields above the scale A. A 
complete survey of the operators generated at dimension five and six can be found in [20]. 

We will be interested only in the leptonic sector, for which the SM Lagrangian is: 

/ 2 \ 2 

4m = LilpL + EilpE - [y 4 LHE+ h.c) +D /l H*D li H- X H (h*H- y J , (2.2) 

1 Which can be ambiguously read as either a fluctuation to set upper limits to the branching ratio, or as 

a signal and a measurement of the same branching ratio. We give more details on section 2.1. 
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where L and E are the lepton weak doublet and singlet fields (triplets in flavor space, 
indexes are suppressed), H is the Higgs doublet, D^ is the appropriate covariant derivative 
and 2/4 is a matrix in flavor space. 

In the EFT, LFV Higgs interactions are mediated at lowest order by the dimension six 
operator 2 * 

£( 6 ) D c ™6 m + h.c = ^LHH^HE + h.c, (2.3) 

A z A z 

where cqh is a matrix in flavor space and we define the operator Oqh = LHH^HE. In 
principle, the couplings for different chirality combinations may differ (e.g. c T ^ R ^ c qj^ r ) 
but, since we will be always looking at initial and final states containing both combinations, 
we consider a single coupling cqhi understood to be the average value cqh = \J 
where LR and RL are the chiralities of heavier and lighter fermions, in this order. This op¬ 
erator, in conjunction with the renormalizable Yukawa couplings ( 2 / 4 ) generates the effective 
flavor violating Yukawas, given in the mass basis by: 

V2H = 2/4 + 4|?. (2.4) 

In our analysis we assume that the diagonal entries of this matrix are close to their SM 
values (i.e., Yu ~ rrii/v). 

In the presence of a non-diagonal Y, a one-loop contribution to radiative lepton decays 
h —•► lj 7 is generated (see Fig. la). If O^h is the only higher dimensional operator relevant 
for Higgs flavor physics, Y,y is constrained by the experimental bounds shown in Table 
1. This diagram however must contain an insertion of the lepton Yukawa coupling, which 
greatly suppresses the radiative decay, allowing the LFV Higgs couplings Y t] to be sizeable 
while still respecting the limits in Table 1. 



Figure 1: (a) One-loop contribution to radiative lepton decay due to off-diagonal Yukawa 
couplings, (b) The LFV decay h —> Tfi. 

Once a large enough is allowed, one may study the Higgs LFV decay h —> lilj 
generated by it (see Fig. lb). In the case of h —> /re, the strong bound of /r —> e 7 (Table 
1) puts this process at a rate beyond the reach of the LHC: BR(h —> /re) < 1CV 8 [3]. The 
authors of [3] claim that the non-diagonal Yukawas for h —>• re and h —> Tfj. are still allowed 
to be big enough to produce a branching ratio at an observable level. Recently, CMS has 
released a data on the direct search for h —> r/r and there is an excess of events of a 2.5 a 

2 There are no 5-dimensional operators not involving the neutrino sector and other dimension six opera¬ 

tors are either irrelevant to flavor violation or equivalent to the operator of eq. (2.3) [3]. 









Process 

Upper bound on BR 

M -► e 7 
r —> eq 

T H 7 

< 5.7 x 10 -13 [21] 

< 3.3 x 10 -8 [22] 

< 4.4 x 10 -8 [22] 

h —> Tfj, 

< 1.57 x 10 " 2 [ 6 ] 


Table 1: Experimental bounds on LFV. 


significance in the channel [6]. That could be read as a statistical fluctuation and be used 
to put a bound of BR(/i —> t/jl) < 1.57% at 95% confidence level, or as a signal, which then 
leads to BR(/i —»• r/r) = (0.891q'| 7 ) % [ 6 ], only with more data will the situation become 
clear. The possibility of confirming this as a signal is specially interesting because, as we 
will see in the following chapters, the conclusion of [3] does not hold if there are other 
sources of flavor violation (besides O^h), and a confirmation would disfavor models with 
such sources. We will focus on this process for the remainder of this text. 

2.2 Interplay with Radiative LFV 

The conclusion of [3], that Higgs LFV decays may be observed in the tau sector, holds if 
the only contribution to the tau radiative decay comes from the Higgs one-loop diagram 
(Fig. la). In general, however, other sources LFV above the scale A may generate the 
dipole operators: 

C EFF D + e^O S7 + h.c, 

= e C j^LHcr a pEF a P + Hcr a pEF a P + h.c, (2.5) 

where F a @ is the electromagnetic field strength, e is the electric charge and the dipole 
operators are defined to O 6 7 = LF[(j a pEF a P and O s 7 = LHH^Ha a pEF a P. Since these 
operators contribute at tree level to the flavor violating radiative decays they cannot be 
immediately dismissed. In the literature, one usually ignores O s 7 on account of it’s larger 
suppression with the scale A [5, 23-25]. However, in cases where Oq^ is not present (i.e., 
it is independently suppressed or unrelated to LFV, see the next section for examples), it 
may become the leading contribution to LFV radiative decays. 

Thus, in general, the bounds from Table 1 will apply to a combination of and 

generally making the restrictions on stronger, unless we restrict the analysis to 
specific UV theories in which and ^ are suppressed in relation to , either by having 
the dipole operators be generated a higher loop orders or by different physics at higher 
scales (making A in eq. (2.5) different and bigger than A in eq. (2.3)). 

It turns out that it is quite hard to find simple UV completions that generate cqh C 6 7 
and cqh cs 7 at a given scale A. 

For the sake of illustration, lets first assume that Oqh is generated by the UV theory 
at loop level (one loop or more). If there is charge flowing through the loop, one can easily 
add an emitted photon, generating C> 8 7 with cs 7 — cqh (see Fig. 2a). The caveat here 
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Figure 2: Relation between the LFV Higgs coupling and higher dimensional dipole op¬ 
erators. (a) Obtaining the dimension eight dipole from the loop generated LFV operator, 
(b) Obtaining the dimension six dipole from the tree generated LFV operator. 


is that one could propose UV theories in which there is no charge going through the loop 
at leading order, in this case cs 7 would be suppressed by loop factors. We provide one 
such example in section 3.4. One could also obtain 0%^ from 0 s 7 by closing a loop with 
two external Higgs lines or by simply removing them. In the first case C6 7 will be a loop 
factor smaller than cs 7 and in the second case the relation between C6 7 and the other two 
couplings depends on the specific couplings of the Higgs with the particles going around 
the loop, making it effectively independent. 

On the other side, if Oqh is present at tree level in the UV theory, one cannot simply 
attach a photon and get Os 7 , because the diagram must vanish by Ward Identities. But 
now, one can obtain Oq^ by closing two of the Higgs legs in a loop and attaching a photon 
to that loop (see Fig. 2b). In this case C6 7 will be smaller than cqh by a loop factor, but 
since it contributes to the radiative decay at tree level while cqh contributes only at loop 
level, the factors cancel and the dipole operator must be taken into account. 

We see that one may connect the dipole operators to Oqh and as such the Wilson 
coefficients C6 7 , cs 7 are in general not independent from cqh- Of course, it is always possible 
to imagine that contributions from additional degrees of freedom to the dipole operators, as 
well as (approximate or exact) flavor symmetries, unrelated to Oqh may greatly suppress 
the values of C6 7 , cs 7 with respect to cqh, allowing the latter to be generated at lower scales 
and to have observable effects. In practice, however, a survey of the recent literature on the 
subject [7-14, 17-19] reveals that in many well motivated UV completions, such as SUSY 
and Composite Higgs scenarios, there is a great correlation between the sizes of ce 7 , cs 7 
and cqh, as we will discuss in the next section. 

If the theory that UV completes the SM is a gauge theory, it is possible to show that the 
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radiative dipole operators are always generated at one-loop level [26]. This follows directly 
from imposing the QED Ward identities in the radiative decay invariant amplitude. Such a 
consideration motivates extracting a universal 1/(167r 2 ) factor from the Wilson coefficients 
C6 7 and cs-y- Here we do not do this, because while these operators are always generated 
at loop level, the Wilson coefficient cqh may be generated at tree or loop level, depending 
on the specific UV theory. As such, we will at first assume that Cqh, cq 7 and cs 7 are 
generated at the same order, implying that any loop factors are the same in both, to obtain 
our bounds. Later, when we discuss specific models, we will note the cases where cqh is 
generated at tree level, or more generally at a lower loop order than the dipole operators, 
and hence is less suppressed by one or more loop factors. 

Assuming that c@ 7 , Cs 7 and cqh me of the same order, we may compare the contribu¬ 
tions to the r —> /ry decay rate generated by each of the three higher dimensional operators. 
In the following expressions we neglect interference between channels, since our goal is to 
find out in what regime only one of the Wilson coefficients is sizeable and the other two are 
small enough to be neglected. The contribution of Oqh is shown in Fig. la and the contri¬ 
butions of Og 7 and O s 7 are obtained by substituting H by in the rightmost diagrams 
of Fig. 2. The expressions are: 


where: 
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( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 


here we take the top Yukawa yx = (vt)sm = tut/v = 0.70. 

We note that in the processes involving lepton radiative decay one must include a 
certain class of two loop diagrams, the so-called Barr-Zee diagrams [3, 27]. These diagrams, 
while of higher loop order, are less chirally suppressed, as such become parametrically more 
important than the diagram in Fig. la and involve the same non-diagonal Yukawa coupling. 
The Barr-Zee contribution is denoted by c 2_loop in eqs. (2.8) and (2.10). We also note that 
we are assuming that the one loop and two loop contributions have a positive relative sign 
i.e., there is no destructive interference between these diagrams. One way a cancellation 
may be achieved is if the top Yukawa (that enters in c 2_loop ) has a negative sign with respect 
to the top mass, yx = —mxV2/v. In this case the numerical value of our bounds is changed, 
but qualitatively our results are unchanged. 

In Fig. 3, we plot the branching ratio BR(r —> juy) calculated using eqs. (2.6) to 
(2.8) as a function of the UV scale A. We see that, for similarly sized Wilson coefficients, 
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O 6 7 dominates completely the amplitude, generating the strongest bound on A, while even 
the highly suppressed dimension eight operator O becomes important for a cutoff below 
approximately 25 TeV and c 6 7 ~ Cs 7 ~ cqh ~ 1- In general Og-y will be less important than 
Oqh for: 

A > (25.7 TeV). F^. (2.11) 

V c 6 H 



A (TeV) 

Figure 3: r —> /ry branching ratio (including the Barr-Zee diagrams) from various opera¬ 
tors as a function of the new physics scale A. The continuous black, dashed blue and dotted 
red lines are respectively the branching ratios generated by © 67 , (^87 an d Oqh- Here, we 
consider ~ cg^ ~ cqh ~ 1 , and interference effects are neglected. 

In the case where Og^ is generated at one higher loop order than Oqh, there is an extra 
factor of ~ 1 /( 167 t 2 ) in the definition of the operator, as discussed before. In this case the 
crossover happens at: 

A ~ (2.0 TeV), (2.12) 

V c 6 H 

It should not be so surprising that a dimension eight operator may play a role in such 
a case, since the diagram of Fig. la that generates the dipole operator at one loop has an 
additional suppression due to the r Yukawa insertion. In fact, it is for this same reason 
that we must include the Barr-Zee two loop diagrams when computing the contribution of 
the operator of eq. (2.3) to the tau decay rate. Hence, Og^ cannot be neglected in many 
cases of interest. 

Once we have estimated the size of the various contributions to the LFV decay process, 
it becomes a relevant question to ask whether there are any interesting models in which 
the connection between the operators generating radiative and Higgs mediated LFV decay 
is broken, allowing the Higgs LFV decay to be seen at the LHC. One such example is a 
type-III two Higgs doublet model, explored in [14], while another one, explored recently in 


- 7 - 





H 

H 

H 

H 

i , H 

T 



M Low _ T 

W M 

Y 

-ft- 

Y 

Y 

Energies 

ic 6H 


Figure 4: Tree-level insertions of the Higgs generating Oqh in the composite sector model 
of [7]. The double line represents vector-like heavy lepton states, while the thin lines are 
chiral elementary fermions. 

[17], is provided by a scalar gauge singlet model, which may be a viable DM candidate. 
We discuss both these scenarios as well as the SUSY and Composite Higgs ones in the 
next section, keeping attention to the relative size of the higher dimentional LFV operators 
produced by each model and whether or not they may generate the h —> Tfi process at an 
observable level for the LHC. 

3 Results for Specific Models 
3.1 Composite Higgs 

In models where the EWSB happens due to the dynamics of a new strong sector, such as 
Randall Sundrum (RS) [28] models and their discrete versions (Quiver or IV-site models [29— 
31]), heavy vector-like leptons may mix linearly with the chiral SM leptons, in the so-called 
partial compositeness scenario. In these models, the Higgs boson is a light composite bound 
state, and does not couple directly (i.e. in the flavor basis) to the SM leptons. Following 
a recent paper on a simplified model containing vector-like leptons [7], one may write the 
Lagrangian 

C = M\iL l ^ r + MX e E R ^> L - Mqff - Mc e ff + Y^ L H^ R + Y^ R H^ L + h.c. (3.1) 

where Lj j and E R are the elementary electroweak lepton doublet and singlet, \L and T 
are respectively the corresponding SU(2)l doublet and singlet vector-like fermions of the 
composite sector, H is the Higgs doublet, M is the mass scale of the heavy resonances, 
A/, X e , Y, Y are matrices in flavor space (all fermion fields are flavor triplets) and we are 
in the {T, \L} basis where q and c e are flavor diagonal. 

In this class of models, Oqh is generated at tree level by integrating out the heavy 
fermions, as per the diagram of Fig. 4. It is essential in obtaining this operator that there 
be non-zero “wrong” chirality couplings Y of the Higgs to the composite sector. Such an 
operator is naturally obtained when considering the brane localized (pure composite) Higgs 
as a limit of a bulk localized Higgs [18]. Upon closing two of the external Higgs legs and 
adding an external photon, one generates and it is thus easy to see that the Wilson 
coefficients cq 7 and cqh are proportional (as in Fig. 2). Furthermore, because the required 
chirality flip occurs on a vector-like fermion line, it does not suffer from chiral suppression by 
the tau Yukawa. Because of this strong correlation, the stringent bounds on the radiative 
decay r —> /ry force the Higgs decay h —> tji to be well below observational level. The 
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authors of [7] obtain BR(/i —> t^l) < 8.6 x 10~ 6 , which is much more restrictive than the 
current bounds obtained by direct searches: BR(/i — > t/i ) < 1.57 X 10 ~ 2 [6]. 

In a more complete model, such as the RS model explored in [32], there are additional 
contributions to and cs 7 , due to composite gauge bosons (electroweak Kaluza-Klein 
(KK) modes). In this case, the main contribution is due to the first charged resonance 
W', which sets the strongest bound on the KK scale, of A ~ 5 TeV for Y, Y of order one. 
Barring cancellations with the typically smaller contributions from the neutral sector (Z' 
and h), one may use such a bound on the KK scale to obtain the expected value of the 
h —> Tfj, branching ratio to be of order m/ l /(87rA 4 ) r\j 10 14 , for c 6H — 1 - 

Both these cases consider the Higgs to be a bound state of the strong sector that is 
accidentally light. A well motivated proposal to explain the lightness of the Higgs with 
respect to the other resonances is that it is a pseudo Nambu-Goldstone Boson (pNGB). 
In this case, one assumes the strong sector has a large global symmetry Q, spontaneously 
broken down to H, such that the NGBs parametrizing the G/'H coset contain the Higgs 
doublet, conveniently represented by the non-linear sigma model held £ = exp (ih/f), 
where h are the components of the doublet H and / is the Higgs decay constant. The 
symmetry Q is then broken by the Yukawa couplings and by gauging a subgroup %’ D 
SU{2)l X U(1)y, giving the Higgs, at the one loop level, a small mass compared to the 
strong sector resonances. 

The requirement that the Higgs couplings must respect a global symmetry of the strong 
sector, in which the SM gauge symmetry is embedded, modifies the effective Lagrangian, 
such that Oqh must be aligned to the dimension four Yukawa coupling, suppressing flavor 
violation [33] 

Following [33], this alignment can be seen by promoting the strong sector operators 
'k, \k to full representations of Q. Then, the Yukawa couplings Aj, A e of eq. (3.1) are made 
formally invariant under Q by embedding the elementary fields in incomplete Q representa¬ 
tions, where the additional degrees of freedom are spurions. In place of the direct couplings 
to H, the strong sector must couple Q invariantly to £. 

Global invariance then forces the couplings to the doublet H to take the schematic 
form 

TJ~ / TJ f TJ \ 

/Aiy ( 2/4 + 2/6 - p— 1 - j Er = L l P{Y)E r , (3.2) 

where P(£) is a polynomial in £, projected over the elementary fields. As long as the 
elementary fields are coupled to a single representation (e.g. the 5 of 50(5) in the MCH5 
model [35], in which case P(£) = ££ T —> sin(/r//) cos(/i//)), then all the higher dimen¬ 
sional non-derivative operators in H are aligned and can be simultaneously diagonalized. 

In this case the leading contribution comes from the kinetic mixing generated by inte¬ 
grating out the heavy leptons, which is typically smaller than (2.3) and leads to even lower 
Higgs LFV decay rates. One may then conclude that the observation of the process h t/j, 
at a rate near the current experimental sensitivity would disfavour this class of models. 

3 While [33] considers the quark sector, the lepton sector is analogous. See [34, 36]. 
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3.2 2HDM 


In models with two electroweak Higgs doublets, the Yukawa Lagrangian for the charged 
leptons is written as 

C = yiL^iE + y 2 L$ 2 E + h.c, (3.3) 

where 4>i and 4>2 are the Higgs doublets and y i, y 2 their corresponding Yukawa matrices. 
Neglecting CP violation for simplicity, the neutral CP even scalar mass eigenstates are 
related to those in the flavor basis by 

h = sin(a — /3)4>i + cos(a — /3)(f> 2 , 

H = cos(a — (3)(J) i — sin(a — (3)(f> 2 . (3-4) 

We consider the case where the lightest eigenstate h is the 125 GeV Higgs. Limits on 
deviations from the SM couplings then demand sin(a — j3) ~ 1, while flavor constraints on 
B meson physics demand m E > 300 GeV [37]. If there is a parity symmetry distinguishing 
the two doublets, (e.g. type II 2HDM), one may forbid the couplings of 4>2 to the charged 
leptons at tree level. This case is closely related to the SUSY models considered in the 
following section and will be discussed there. 

In the general case where there is no such parity, it is generally impossible to diagonalize 
simultaneously y\, y 2 and the Higgs mass matrix, and we then see that a Higgs LFV 
operator is already present in the renormalizable Lagrangian in the Higgs mass basis. In 
order to connect this result with our effective theory framework, one must integrate out the 
other physical scalars (assumed heavier than the 125 GeV Higgs). Since, by eq. (3.4) the 
physical scalar h lives in both doublets, the contributions coming from 4*2 proportional to 
y 2 cos(a — (5 ) to the Yukawa matrix do not decouple for big m E and are not suppressed by 
the new physics scale, being in general stronger than the higher dimensional operators. 

Recently, the authors of [14] have shown that in such a case one may obtain the h —>• ry 
process at an observable level, despite the limits from radiative decays. The reason being 
that the Higgs LFV is generated by a renormalizable operator, at tree level, while the tau 
radiative decay is loop suppressed. Furthermore, in certain regions of parameter space, 
a cancellation between the heavy scalar and pseudoscalar becomes possible, leading to 
additional suppression of the radiative decay. 

3.3 SUSY 

In supersymmetric models such as the Minimal Supersymmetric Standard Model (MSSM), 
anomaly cancellation and holomorphy force one to include two Higgs doublets of oppo¬ 
site Hypercharge, that give mass separately to up type and down type fermions. This 
corresponds then to a type II 2HDM, at least at tree level. 

At one-loop order, one may generate Yukawa couplings with the “wrong” Higgs boson, 
as in the diagrams of Fig. 5. Consider the soft SUSY breaking trilinear terms, and slepton 
masses 

4oft = ~A e ELH d + h.c. - M 2 l \L\ 2 - M 2 e \E\ 2 , (3.5) 

where E, L are the left and right slepton fields and H d is the down type Higgs doublet. 
The trilinear couplings A e and the soft masses Ml, M e are matrices in flavor space that 
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Figure 5: Diagrams leading to coupling between charged leptons and the “wrong” Higgs 
doublet (i.e., the one not responsible for their mass at tree level). The slepton, higgsino, 
wino and bino fields are denoted l , H 'j*’ 2 , and B° , respectively, and the cross denotes 

a mass insertion. Chirality labels are suppressed for simplictity. 


in general are not aligned both with each other and the Yukawa matrices. This allows the 
diagrams of Fig. 5 to generate small flavor violating transitions. Adding a photon to the 
loop, and exchanging the Higgs for its vev, one may then generate the radiative lepton 
decay with a Wilson coefficient that is directly proportional to the Higgs LFV one. Just as 
in the composite Higgs case, both Wilson coefficients, cqh and C6 7 are correlated, however, 
unlike in that scenario, they are both generated at the same loop order 41 . 

Since O 6 7 contributes to r —> pry at tree level, C6 7 will strictly bounded. The relation 
between cqh and C6 7 then implies BR (h —> r/x) to be very small. The authors of [10], for 
instance, obtain BR(/i —> t/j.) < 10 -4 from BR(r — > /xy) < 3.1 x 1CV 7 (which is an older 
constraint, superseded by the one shown in Table 1). A R-Parity violating scenario was 
explored in [11] with BR(/i —y r/x) < 1CV 4 5 as a constraint. 

This puts the LFV Higgs decays at least two orders of magnitude below the current 
experimental sensitivity, and disfavors this class of models in case the current excess turns 
out to be signal. 

3.4 Higgs Portal 

One class of models where radiative lepton decay may be small enough to satisfy bounds 
and still allow observable Higgs LFV decays is given by the Higgs portal models, where a 
scalar singlet odd under a Z 2 symmetry is added to the SM. The coupling terms of this field 
to the Higgs and leptons, as well as its self couplings, are given by the effective Lagrangian: 

4 While the diagrams of Fig. 5 are directly related to the LFV Higgs decays, there are other classes of 

diagrams that contribute to r —> /xy, as in [10], and many more if R-Parity is violated, see [11]. 
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Figure 6: LFV diagrams for a Higgs Portal model, (a) Higgs LFV interaction at one loop, 
(b) Radiative lepton flavor violation at two loops. 


C S H = \s 2 H'H - ^S 2 LHE - ^S 2 + (3.6) 

where S is the singlet field, £, y are coupling constants ( y is a matrix in flavor space), and 
A is a scale associated with some heavier sector (e.g., vector-like leptons or an additional 
scalar weak doublet) which generates the effective interaction of S with SM fermions, not 
to be confused with the mass of S, ms • The flavor violating decays are generated by the y 
couplings, which may be misaligned with respect to the SM Yukawas. 

The mass of S is given at tree level by 

m S = m S 0 + 2 V ' ( 3 - 7 ) 

At one loop order, there are additional contributions from A s and y which we neglect. 

In this class of models, the Higgs LFV interaction is generated at one loop order by 
the diagram of Fig. 6a, which becomes Oqh if we also integrate out S. 

The radiative lepton flavor violation is generated only at two-loop order by the sunset 
diagram of Fig. 6b. Both and Os-y are obtained in the large ms limit, O s 7 is directly 
obtained from Fig. 6b, while Og 7 is obtained by closing a third loop with two of the three 
Higgs lines in that diagram, making this operator even more suppressed. 

We note that, while from a low energy point of view, ms is the cutoff scale that 
determines the suppression of the irrelevant operators, an explicit calculation from the 
diagrams of Figs. 6a and 6b shows that cqh an d cS'y depend on ms only logarithmically, 
becoming important for a larger gap between ms and the scale A, see below ’. 

Because the radiative decays are generated at two loops, they are naturally suppressed 
with respect to the Higgs LFV rates. Furthermore, since S is uncharged, it is impossible to 
add a photon to the loop of Fig. 6a, and the Wilson coefficients cqh and Cs 7 are uncorrelated. 
In this case, we find that the conclusions of [3] apply. 

5 For a large gap, of order ms /A < 0.37, perturbativity is lost and one must resum the large logarithms 
with the RG. This is beyond the scope of our present analysis. 
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A=410 GeV and y r „=47T 


A=410 GeV and BR(h->r^)~10' 2 



Figure 7: Higgs portal coupling £ as a function of the scalar mass ms- We show contour 
lines for (a) BR(/i —> r/r) with y T ^ = 47T and for (b) y Tfl with BR(/i —> r/r) = 10 -2 . The 
region A = 410 GeV and y Tfl < 47r satisfies perturbativity bounds (see text). 


We estimate the size of the contributions to h —)• ry and r —» fi'y to be, respectively: 


c <Wj2 = 


^87 


ev 


A 4 


gj/r/j V 2 
32tt 2 A 2 




— 1 sin 1 


yr^ykkm T v 2 2 m\ 

6 2(16vr 2 ) 2 A 4 g A 2 ’ 



(3.8) 


where the result for cs 7 is only a leading order estimate. 

Using eqs. (3.8), we can put bounds on the parameter space of this simple model. The 
constraints are that the radiative r decay be below the experimental rate of Table 1, that 
the gap between ms and A, as well as the size of the couplings £ and y Tfl be small enough 
to mantain perturbative control. Explicitly, we use the perturbativity contraints 

, / 2 A 2 A 2 

;$ 327T 2 —2 ~ 315.8— 2"5 

1 > X ~ E ~ l ~ °' 37 ' ( 3 ' 9 ^ 

In the second equation, E = 2.718 ... is the Euler number. We note that no stringent bound 
on y Tfl is implied by cs 7 , because of the r Yukawa suppression. Imposing these constraints, 
we obtain the plots of Fig. 7. 

One sees that, satisfying the constraints, it is possible to generate a Higgs LFV decay 
rate of about BR(h —> T/i) = 10~ 2 , observable at the LHC, while the radiative tau decay 
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is well under control, with a rate of order BR(t fj.'j) < 1.5 x 10 -1 ', where the maximum 
value is obtained for couplings saturating the perturbativity bound. However, if S was 
stable (being the lightest Z 2 odd particle), it’s thermal relic abundance would have to be 
below or equal to the observed Dark Matter abundance, fig < £Idm = 0.227, implying 
£ < 3 x 10 _ 4 (ms/l GeV) [38, 39], which is far too small to allow cqh to be sizeable. One 
way to avoid this is to extend the scalar sector to obey a larger global symmetry than Z 2 . 
In this case, only the lightest scalar must satisfy the DM constraints, while it’s heavier 
partners may have larger couplings to the Higgs, of order £ ~ 10 and so allow the h — > t^i 
rate to be generated at an observable level. 

This class of models then allows one to evade the bound on r —> fi'y and still generate 
h T[i at a rate near the current experimental sensitivity. A signal of this decay would 
then be indicative that this kind of model may be playing a role. 

Detailed phenomenology of the scalar sector and the viability of the lightest scalar as 
a DM candidate will be presented elsewhere, however, see [17]. 

4 Conclusions 

We have examined the possibility of a new physics contribution to the Higgs LFV process 
h —> T/JL, which has been recently constrained by CMS to BR(h —> rfj,) < 1.57% at 95% or 
that may be rather seen as a signal, with BR(h —> T{i) = (0.89^0 3 ?) % [ 6 ]. 

From a model independent perspective, in order to obtain a signal big enough to be 
detected or restricted by the LHC (BR(h — > t/i ) ~ ICR 2 ) one must consider models in which 
the Higgs LFV operator Oqh is sizeable but does not violate the important constraint from 
radiative lepton decay, BR(r —> ^ 7 ) < 4.4x 10 -8 [22], This usually requires the contribution 
of the dipole operators O 6 7 , Os 7 to the tau decay rate to be negligible, which implies that 
they are generated at a higher loop order than Oqh by the UV continuation of the SM. 
The dipole operators, specially Os 7 , have been generally neglected in model independent 
approaches to LFV Higgs decays. We show that the contributions of these operators can 
be important and thus put strong bounds on Oqh- 

Looking at specific classes of models found in the literature, we see that typically the 
Wilson coefficients of CVp Os 7 are correlated with Oqh , and make it difficult to avoid the 
radiative decay bound and still get Higgs LFV at an observable level. This is the case of 
SUSY, where Oqh is generated at one-loop order, and Composite Higgs models, in which 
it is obtained at tree level. In both these cases, one may add a photon to Oqh and thus 
get with a similarly sized Wilson coefficient. For this reason, these kinds of models are 
disfavored by the data, should the measured rate for h —>■ r/i be confirmed as a signal. 

On the other hand, there are certain kinds of models in which one may avoid this con¬ 
clusion. Examples are a type III Two Higgs Doublet Model, where Higgs LFV is produced 
by a renormalizable operator, and models with an extended, gauge singlet scalar sector, in 
which case the tau radiative decay only appears via Os^, generated at two loops, and is 
negligible. These options for extending the Standard Model may be favored, should this 
signal persist. 
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More generally, any model aiming to explain the data should have a mechanism to 

suppress r —> /ry independently from Higgs LFV, by guaranteeing that the dipole operators 

are generated at a higher loop order than Oqh- 
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